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ABSTRACT 

Accretion and star formation (SF) rates in low redshift SDSS type-II active galactic 
nuclei (AGN) are critically evaluated. Comparison with photoionization models indi- 
cates that bolometric luminosity (Lboi) estimates based on L([0 in] A5007) severely 
underestimate L uo i in low ionization sources such as LINERs. An alternative method 
based on L(H/3) is less sensitive to ionization level and a novel method, based on a com- 
bination of L([0 in] A5007) and L([0 1] A6300), is perhaps the best. Using this method 
I show that low ionization AGN are accreting faster than assumed until now. Signifi- 
cant related other findings are: 1. Any type-II AGN property related to the black hole 
(BH) mass is more reliably obtained by removing blue galaxies from the sample. 2. 
Scyfcrt 2s and LINER 2s form a continuous sequence of L/L Edd with no indication for 
a change in accretion mechanism, or mode of mass supply. There are very few, if any, 
LINERs in all type-I samples which results in a much narrower L/LEdd distribution 
compared with type-II samples. 3. There is a strong correlation between SF luminos- 
ity, Lsf, and L uo i over more than five orders of magnitude in luminosity. This leads to 
a simple relationship between bulge and BH growth rates, g(bulge) / 'g(BH) oc L^' 2 , 
where g{bulge)/g{BH) ~ 115 for L fcoZ =10 42 ergs s _1 . Seyfert 2s and LINER 2s follow 
the same LsF-Lboi correlation for all sources with a stellar age indicator, D n 4000, 
smaller than 1.8. This suggests that a similar fraction of SF gas finds its way to the 
center in all AGN. 4. L uo i, Lsf, L/LEdd and the specific SF rate follow D n 4000 in a 
similar way. 
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1 INTRODUCTION 

Black hole (BH) masses in thousands of type-I active galactic 
nuclei (AGN) can now be obtained from optical-UV spec- 
troscopy. The method is based on the known relationship 
between the luminosity of the non-stellar continuum at some 
wavelength (e.g. 5100A) and the mean broad line region 
(BLR) size derived from reverberation mapping (RM; e.g. 
Kaspi et al. 2000; Kaspi et al. 2005; Vestergaard and Peter- 
son 2006; Bentz et al. 2009). This size is combined with a 
mean gas velocity estimated from broad emission line pro- 
files (e.g. H/3), and the virial assumption about the gas mo- 
tion, to obtain Mbh. The accuracy of the method has been 
discussed in various papers (e.g. Bentz and Peterson 2006) 
and is estimated to be a factor of about 2. Complications due 
to the effect of radiation pressure force on the motion of the 
BLR gas (Marconi et al. 2008) are probably not very impor- 
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tant, at least for low redshift low luminosity AGN (Netzer 
2009, hereafter N09; see also Marconi et al. 2009). Estimates 
of the normalized accretion rate, L/LEdd, are obtained by 
combining the derived Mbh with estimates of the bolometric 
luminosity, Lboi- 

Different methods are required to obtain mass and ac- 
cretion rates in type-II AGN, where the non-stellar contin- 
uum is not directly observed. In low redshift type-II sources 
situated in bulge dominated hosts, the stellar velocity dis- 
persion in the bulge of the galaxy, a,, can be transformed 
to BH mass (e.g. Tremaine et al. 2002). The bolometric lu- 
minosity is usually estimated from known relationships be- 
tween Lboi and the luminosity of certain narrow emission 
lines. The assumption is that the line luminosity represents 
the same fraction of L^oi in all AGN. A line that is com- 
monly used is [O in] A5007 and the scaling of its luminosity 
to L 00 i is obtained from observations of type-I sources, where 
both the line and the non-stellar continuum are directly ob- 
served. Such scaling is discussed and explained in several 
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papers, e.g. Heckman et al. (2004), Netzer et al. (2006), 
Kewley et al. (2006; hereafter K06), Kauffmann and Heck- 
man (2009; hereafter KH09) and N09. A related method 
which is less sensitive to reddening correction is based on 
the luminosity of the mid-lR line [O/V] 25.9/im (see discus- 
sion and suggested conversion factors in Dasyra et al. 2008) . 
Unfortunately, the number of sources with such mid-IR mea- 
surements is very small. Yet another method is based on the 
luminosity of the hard X-ray continuum which is directly ob- 
served in most type-I and type-II AGN. While potentially 
insensitive to reddening and other complications, the exact 
bolometric correction factor required to convert the 2-10 
keV luminosity to Ltoi depends on the global SED and is 
still debatable (see e.g. the differences between Marconi et 
al. 2004 and Vasudevan and Fabian 2007). 

The emission line scaling method was used by Heckman 
et al. (2004) to investigate mass and accretion rates in a large 
sample of data release one (DR1) Sloan digital sky survey 
(SDSS; York et al. 2000) type-II AGN. L([0 in] A5007) used 
in this case included no correction for reddening and the 
conversion assumed Ltoi — 3500L([O in] A5007). The more 
recent KH09 analysis is based on reddening corrected values 
of L([Oiii] A5007). This method reduces the scatter due to 
reddening and accounts, more reliably, for known differences 
in reddening between galaxies of different stellar popula- 
tions (K03, K06). The reddening-corrected scaling adopted 
by KH09 for all type-II sources is L bo i= 600L([O ill] A5007). 

This paper addresses BH accretion rates, star formation 
(SF) and SF rates (SFR) in low redshift type-II AGN. The 
aims are to test various claims about the accretion mecha- 
nism in LINERs, the dependence of BH accretion on host 
galaxy properties, and the correlation between Ltoi and SF 
luminosity (Lsf)- This requires re-evaluation of the meth- 
ods used to estimate Ltoi in type-II AGN which is presented 
in §2. A new, improved Ltoi indicator is used, in §3, to re- 
evaluate Mbh and L/Lem in LINERs and in Seyfert 2s. 
I also use the same indicator to test various AGN and SF 
correlations. §4 summarizes the main results of the paper. 
Throughout the paper I assume standard cosmology with 
H =70 km/sec/Mpc, fi m = 0.3 and fi A = 0.7. 



2 L([0 in] A5007) L([0 i] A6300) AND L(H/3) AS 
BOLOMETRIC LUMINOSITY INDICATORS 
IN TYPE-II AGN 

2.1 Theoretical line and continuum luminosity 
ratios 

The current analysis applies to narrow emission lines in 
galaxies that contain both an active AGN with a narrow 
emission line region (NLR), and starburst (SB) ionized gas. 
The term SF is perhaps more appropriate to discuss the var- 
ious processes addressed below but SB has been used in the 
past to describe such galaxies and I keep to this terminol- 
ogy when discussing those issues where the name SB is com- 
monly used. While the NLR properties have been studied, 
extensively, observationally and theoretically, the combina- 
tion of AGN and SB excited gas in the same host presents 
a real challenge, especially in sources observed with a large 
entrance aperture. A main tool for distinguishing the two 
types of excitation is based on BPT (Baldwin, Phillip and 



Terlevich 1981) line ratio diagrams where AGN and SB re- 
gions are well separated. Detailed explanation, and critical 
discussion of such methods can be found in numerous papers 
including Veilleux and Osterbrock and (1987), Kauffmann 
et al. (2003a, hereafter K03), K06, Brinchmann et al. (2004, 
hereafter B04), and Groves et al. (2006a; 2006b). 

Photoionization by a non-stellar radiation field can ex- 
plain the large observed range of ionization and excita- 
tion conditions in AGN. In particular, it can explain the 
L([Om] A5007)/L(H/3) (hereafter [OlIl]/H 3 ) line ratio in the 
highest ionization type-II Seyfert galaxies ([Olll]/H^~ 10) 
as well as in type-II LINERs ([O iii]/H,3~ 1). The differences 
are attributed to changes in the ionization parameter of the 
line emitting gas (e.g. Ferland and Netzer 1983; Netzer 1990; 
Groves et al 2004). Models invoking shock excitation have 
also been considered (see reviews in Maoz 2007; Ho 2008). 
Here I consider photoionization as the only ionization and 
excitation mechanism for the AGN gas. I use S2 to refer to 
high ionization Seyfert 2s, L2 to describe type-II LINERs 
and AGN to describe the two classes together. 

Assuming ionization by the AGN continuum, it is 
straightforward to calculate the narrow emission line spec- 
trum using standard photoionization models. Realistic mod- 
els must take into account the expected range of density, 
column density, composition, dust content, the spectral en- 
ergy distribution (SED) of the central source and the spatial 
distribution of the gas. Examples are the locally optimally 
emitting clouds (LOC) model of Ferguson et al. (1997), the 
multi-cloud model of Kaspi and Netzer (1999) and the con- 
stant pressure model of Groves et al. (2004). However, the 
calculations of L([0 in] A5007), L(H/3) and L([0 1] A6300) are 
relatively simple and single cloud models represent well the 
mean conditions in most NLR models. This has been il- 
lustrated in several recent papers such as Baskin and Laor 
(2005) and Melendez et al. (2008). 

New single cloud models have been calculated to in- 
vestigate the dependence of the narrow line luminosities on 
various parameters, in particular the continuum SED, the 
ionization parameter, the gas composition and the dust con- 
tent. The calculations were performed with the code ION 
which is a standard, well tested photoionization code suit- 
able for calculating NLR and BLR AGN models (Kaspi and 
Netzer 1999; Netzer 2006 and references therein). Regard- 
ing SEDs, two different possibilities ought to be considered. 
The first is a "standard" SED that represents high ioniza- 
tion AGN and assumed to retain its shape in phases of low 
and high accretion rates. Given the unobserved far-UV con- 
tinuum, the range of possible shapes is large and various 
components must be considered. These can be described by: 
1. The shape of the 0.1-1/xm continuum. 2. The bolometric 
correction factor, BC, which specifies the ratio of the total 
luminosity to the optical continuum luminosity. The BC is 
not directly observed and its value is estimated from anal- 
ysis of various emission lines. 3. The normalization of the 
optical and X-ray fluxes, e.g. the value of a ox . 4. The shape 
of the X-ray continuum. 

Experimenting with a range of possible continuum 
shapes shows that the two that are shown in the left panel of 
Fig.[T]bracket most realistic SEDs. The first of these (SED-1) 
is a combination of a disc-like optical-UV continuum, com- 
bined with a L v oc v~ - 9 X-ray (0.1 < E < 50 keV) power- 
law. In this case, BC, which is defined here as Ltoi/isioo, 
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Figure 1. Left panel: Two assumed SEDs representing the range 
of continuum shapes in AGN. The top curve (SED-1) is made of 
a disk-like optical -UV continuum and a powerlaw F v oc u~ ' 9 
X-ray (E> 0.1 keV) continuum. For this SED, BC=7 and a ox = 
1.35. SED-2 has the same optical and X-ray shapes and the two 
are combined with a F v oc z/ -1 - 5 powerlaw. Here BC=14 and 
olox = 1.05. Right panel: High accretion rate SED (VI) and 
four assumed low accretion rate SEDs obtained by truncating 
the "blue bump" in different ways keeping a ox unchanged (1.05). 



where L5100 is \L\ at 5100A, is 7 and a ox = 1.35. Such 
SEDs are typical of high luminosity AGN. The second case 
(SED-2) combines the above X-ray powerlaw and the opti- 
cal disc-like continuum with a Lyman continuum powerlaw 
given by L„ oc v~ x ' h . In this case BC=14 and a ox = 1.05. 
Such values of BC and a ox are more typical of low luminos- 
ity AGN. The two continua shown in Fig. [T] are normalized 
to have the same ionizing luminosity to emphasize the fact 
that SED-1 produces relatively more photons close to the 
Lyman edge. 

The second possibility is that different phases of activity 
are associated with different SEDs. For example, a luminous, 
fast accreting phase may be associated with an SED similar 
to one of the continua shown in the left panel of Fig. [T] while 
a less luminous, slower accreting rate phase, might show a 
different SED. This may be the result of a change in the 
spectral properties of the central accretion disc, in particu- 
lar a significant reduction in the "blue bump" flux. Such a 
possibility is illustrated by the various curves shown in the 
right panel of Fig. [1] The top SED, marked VI, represents 
the high accretion rate phase. It is similar to SED-1 except 
that a relatively stronger X-ray continuum, typical of low 
luminosity AGN, has been used. In this case BC=12.4 and 
oiox = 1-1- The additional curves, marked V2-V5, repre- 
sent several possibilities where the accretion disc continuum 
is significantly reduced. Such spectra have been speculated 
for LINERs and the ones shown here were motivated by 
the recent work of Maoz (2007). According to Maoz (2007), 
there is a strong UV (2500A) point -source continuum in sev- 
eral nearby LINERs. In those cases, a ox ~ 1, very similar 
to the values observed in Seyfert 1 galaxies. It thus seems 
that either there are no significant changes in the near UV- 
SED between low and high accretion rate phases, or else 
such changes only affect the part of the continuum between 



~ 2500 A and ~ 1 keV). The various curves in the diagram 
illustrate several such possibilities. The SEDs are truncated 
at energies from 0.5 Rydberg, just beyond the 2500A mea- 
sured point, to 1.81 Rydberg, the Hel ionization edge. The 
transitions from SED-V1 to any of the other SEDs mimic 
several possible changes that are consistent with both the 
Maoz (2007) observations and the assumption of a reduced 
disc emission. 

The constant density models considered here assume 
column density of 10 21 ' 5 cm -2 , enough to make the cloud 
optically thick to the ionizing continuum radiation, and hy- 
drogen number density of nn = 10 3 cm" 3 , low enough to 
avoid strong collisional de-excitation of the [O ill] A5007 line. 
These are quite standard in NLR models and are in agree- 
ment with many observations (e.g. Baskin and Laor 2005; 
Melendez et al. 2008; Ho et al. 1997). I calculated various 
line ratios, for all SEDs, over a large range of the hydrogen 
ionization parameter, U = Q(H)/inr 2 nHC, where r is the 
cloud central source distance and Q(H) the rate of emission 
of Lyman continuum photons. This was done for two generic 
types of models, one with ISM-type depletion and one for 
dust- free gas. The gas composition was changed between one 
and four times solar and the depletion fraction is assumed to 
be independent of abundance. The models shown here are 
all for dusty clouds which are thought to be more appro- 
priate for NLR conditions. I have included several models 
similar to the cases discussed in Groves et al. (2004; 2006b) 
where radiation pressure force operating mostly on grains, 
determines the internal pressure in the cloud. Experimenting 
with the various parameters confirms that U, the dust con- 
tent and the SED shape are the most important parameters 
that determine the interesting line ratios. 

The results of some of the calculations are given in 
Fig. The figure shows the fraction of Ltoi emitted by 
the three lines in question, [O 111] A5007, [O 1] A6300 and H/3. 
The ratios are calculated assuming a covering factor of unity 
(fl = 47r). The left panel of Fig. [5] shows a noticeable dif- 
ference between the behaviors of [O in] A5007 and H/3. The 
intensity of [O ill] A5007 depends strongly on the level of 
ionization and L([Om])/Lboi varies by about 1.7 dex over 
the relevant range of 10" 4 < U < 3 x 10~ 3 . In contrast, 
L(H^)/Lboi is basically constant (less than 0.3 dex) over the 
same range which reflects the constant fraction of the Ly- 
man continuum photons absorbed by the optically thick gas. 
Such results are well known and well documented in numer- 
ous publications. In dust-free large U models (not shown 
here), L(H^)/Lboi and L([0 lll])/Lboi are larger by about 
a factor two compared with the dusty models shown here. 
Dusty gas clouds produce less [O ill] A5007 and H/3 photon 
due to the absorption of part of the ionizing radiation by the 
dust and the destruction of locally emitted line photons by 
grains. The first process is more important in high ionization 
gas. The latter affects H/3 more than [O ill] A5007 in lower 
ionization cases since, in this case, the [O ill] A5007 line is 
formed near the illuminated face where escape out is easy. 

Comparison of the models with different SEDs in Fig. [2] 
shows little difference in L([0 lll])/Lboi since the spectral 
shapes of the two are quite similar at high energies. This is 
not the case for L(H,a)/Lboi because of the larger fraction of 
ionizing photons just beyond the Lyman edge in SED-1 (see 
Fig. 1). All this, again, is well known from earlier calculations 
(e.g. Netzer and Laor 1993). 
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Figure 2. Left panel: Calculated L([0 in])/L bol (black) 
L([Ol])/Lboi (blue) and L(H^)/Lboi (red) ratios for various pho- 
toionized, dusty solar composition clouds with full covering of the 
radiation source. Thin lines: SED-1. Thick lines: SED-2. Thick 
dashed line: SED-2 and radiation pressure dominated clouds. 
Right panel: Similar calculations for the VI continuum. Same 
colour coding, same ratios as left panel but vs. [Oi]/[Oiii]. Thin 
lines, solar abundances. Intermediate width lines, 3xsolar abun- 
dances. Thick lines, 3xsolar, constant pressure model. 



I also calculated L([Ol])/Lboi curves and show them 
on the same scale. This ratio is somewhat more sensi- 
tive to U than L(H ( 3)/Lboi but is much less sensitive 
than L([Olll])/Lboi- The differences between [Ol] A6300 and 
[O in] A5007 are clearly seen in the right panel where I plot 
the same ratios as a function of [Ol]/[Om]. This line ratio 
varies with U in a way which is different from [O Illj/H^. All 
this is the basis for the new calibration method discussed in 
§2.2. 

In conclusion, the fraction of Lboi re-radiated by the 
narrow [O ill] A5007 line depends, strongly, on the level of 
ionization of the gas. The fractions of Lboi re-radiated by H/3 
and [O i] A6300 are almost independent of U . L(H/3) seems to 
be the most reliable Lboi indicator for pure AGN excitation. 

Similar line-to-continuum luminosity ratios were also 
calculated for the scenario illustrated in the right panel of 
Fig. [T] Table [1] lists some of the results. The first row in 
the table lists luminosity ratios for a high S2 accretion rate 
phase assuming the VI SED. I fixed the value of U to give 
[O iii]/H,3= 10, typical of many such sources. The other four 
rows represent low accretion rate phases with the V2-V5 
SEDs. Given an SED, I fixed U to produce [Oll^/H^ 1, 
typical of many L2s. 

Examination of the table confirms that changes in SED, 
even as extreme as the ones considered here, do not change 
much the previous conclusions. The line-to-continuum ra- 
tios change in a way similar to the one shown in Fig. [2] such 
that L([0 lll])/Lboi drops much more (factors of 10-20) com- 
pared with the drop in L([0 l])/Lboi and L(H / 3) /Lboi (factors 
1.5-2). The reason is that the difference in ionization poten- 
tial between hydrogen and +1 is small and changes in the 
flux of the +1 ionizing photons are associated with sim- 
ilar changes in the number of Lyman continuum photons. 
Most of the reduction in [Oiii]/H,3 is due to changes in U 



and cannot be explained by changes in continuum shape, at 
least not for the SEDs considered here. 

To summarize, large changes in [Oiii]/H,3 are associ- 
ated with large changes in L([0 lll])/Lboi- In contrast, the 
changes in L(H / 3)/Lboi and L([Ol])/Lboi following a compa- 
rable change in [OiiiJ/H^ are much smaller even for dra- 
matic changes in SED shapes like the ones considered here. 

2.2 Ltoi estimators 

2.2.1 AGN samples BH mass estimates and reddening 
corrected line intensities 

The most reliable estimates of L([0 m])/Lb i and 
L(H,3)/Lboi are obtained from observations of type-I 
AGNs. In such sources, L5100, L([0 in] A5007) and, to a 
lesser extent L(H/3) (due to blending with the broad H/9), 
are directly observed and luminosity dependent bolometric 
corrections are known with reasonable accuracy (Marconi 
et al. 2004; Netzer et al. 2007; Shen et al. 2008; Vestergaard 
2008; see summary in N09). This allows a direct conversion 
of line to continuum luminosity. 

To quantify this, I used type-I data from the SDSS/DR5 
sample described in Netzer and Trakhtenbrot (2007, here- 
after NT07). As explained in N09, the sample is incomplete 
at 2 ^ 0.1 due to the selection method of type-I AGN in 
SDSS data. The completeness level, assuming a flux limited 
sample, is close to 100% at 2 ^ 0.2. The sample contains 
about 9000 radio quite (RQ) type-I AGN with 2 ^ 0.75. 
The NT07 work resulted in line and continuum measure- 
ments, including L([Om] A5007), L(H/3) and L5100, for most 
of these objects. The number of sources in the redshift in- 
terval 0.1-0.2 is 1333 and the number of sources with re- 
liable narrow [Om]/H/3 measurements is 1172. The fitting 
procedure used to measure the lines, and the other details 
of the analysis are described in NT07. The bolometric cor- 
rection factor I used is somewhat different than the one 
used by NT07 and is taken to be BC=9 — logZ/44, where 
La =L 5 ioo/10 44 ergs- 1 (see N09). 

Several type-II samples are used in the present analy- 
sis. The first was extracted from the SDSS/DR4 (Adelman- 
McCarthy et al., 2006) archive which is available on the 
MPA URL siteQ- The archive includes redshifts, stellar ve- 
locity dispersion measured through the 3" SDSS fibre, line 
fluxes, the D n (4000) index (luminosity weighted mean stel- 
lar age) and various other properties. The stellar velocity 
dispersion was used to derive BH mass. This is a reliable 
procedure for elliptical and bulge galaxies but gives poorer 
results for disk and pseudobulge systems. In this work I am 
studying mass and accretion rate distributions in type-II 
AGN and I also make a comparison with type-I sources. For 
the first part I will try to avoid galaxies with less reliable BH 
mass estimates. This requires the removal of blue galaxies 
from the sample and is explained in §3.1. For the compar- 
ison with type-I AGN, I retain the entire sample since no 
separation of blue and red galaxies has been applied to the 
type-I SDSS sample used here. 

The line fluxes for the type-II sources were re-extracted 
from the raw data in the archive and reddening corrected 

1 www.mpa-garching.mpg.de/SDSS/DR4/ 
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Table 1. L\> i estimators for full covering dusty NLRs with solar composition and various SEDs 
shown in Fig. [T] 



SED 


\ogU 


[OmJ/Hfl 


L 6oi /L([Om]) 


L 6oi /L(H/3) 


L ftoi /L([0 i] A6300) 


BC 


VI 


-2.58 


10 


74 


776 


1047 


12.4 


V2 


-3.56 


1.06 


851 


891 


933 


11 


V3 


-3.57 


1.1 


912 


1000 


525 


10.1 


V4 


-3.49 


1.04 


1585 


1660 


501 


9.1 


V5 


-3.33 


1.1 


617 


692 


708 


11.6 



fluxes were obtained in two ways, using the observed H^/H^ 
line ratios. The first assumes standard galactic reddening 
curve combined with the assumption that the intrinsic ratio 
is H a /H/3=2.86. This value is somewhat smaller than the one 
expected for S2s (3-3.1) but is a very good approximation 
for L2s. The second is aimed to duplicate the values used 
in K03, B04 and in Groves et al. (2006a, 2006b). Here the 
extinction is given by A\ cx A~ 0,7 . Part of the motivation to 
use such extinction is related to the issue of foreground vs. 
internal dust. Here it is only taken as one of two possible red- 
dening laws. The A -0 ' 7 reddening law is very flat compared 
with the galactic law and results with significantly larger Ay 
and intrinsic line luminosities. For example, the median red- 
dening for the high S/N (see below) sample of L2s and S2s, 
assuming galactic reddening, is Av = 1.07 and the equiva- 
lent number for the A -0 ' 7 law is Ay — 1.80. Fortunately, a 
consistent use of the two reddening laws (see §2.2.2) results 
in only a small difference in the deduced L^, i. 

Two additional complications related to reddening must 
be considered. First, the emission line reddening in type-II 
AGN is thought to be larger than in type-I AGN. This has 
been discussed in various papers, e.g. Netzer et al. (2006), 
and Melendez et al. (2008). Any reddening correction based 
on average values obtained from type-I samples must in- 
clude this effect. Second, the average Ha/H/3 in previous 
optically selected and X-ray selected samples (e.g. the Bas- 
sani et al 1999 sample) is larger than the one measured here. 
This is possibly explained by K03 who noted the large range 
in Ha/H/3 and the tendency to detect the highest values 
in galaxies containing more SF gas (e.g. smaller D n 4000 
systems). Several of the pre-SDSS selected samples might 
preferentially pick galaxies of this type. The present work 
uses reddening corrected values based on the mean observed 
EW([Om] A5007) in type-I SDSS sources. It is therefore ap- 
propriate to use the mean value of Hcv/H/3 in the SDSS S2 
sub-sample in order to choose the correct scaling between 
type-I and type-II sources. 

The total number of L2s and S2s in the sample con- 
sidered here is about 85,000 and the maximum redshift is 
about 0.25. The main purpose is to obtain the most reli- 
able Lboi by using various emission lines. Therefore, most of 
the results pertain to sources with S/N^ 3 in the Ha, H/3, 
[Oni]A5007, [Nil] A6584 and [Oi] A6300 lines. This reduces 
the sample size to about 42,000. Out of these, about 10,000 
are classified as blue galaxies (§3.1) and are removed from 
part of the analysis, as described below. 

Three sub-samples were created. The first is defined by 
the [Nll]/H Q vs. [Olll]/H,3 K03 criteria for separating AGN 
from SB galaxies. This group is referred to here as the "Ka03 
sample". The second and third groups are obtained by us- 



ing the Ke06 division lines based on the [Ol]/H a and the 
[Nll]/H Q vs. [Oiii]/H/3 line ratios (see K06). These crite- 
ria are aimed to isolate, as much as possible, pure AGN 
from starburst and composite sources. The objects that are 
thought to best represent this group are in the region of the 
BPT diagram defined by the combination of the two K06 cri- 
teria. There are 11,803 such sources and they are referred to 
here as the "Ke06 sample". The S2 sub-group contains 6641 
sources and the L2 group 5162 sources. Obviously, the ratio 
of the populations depends on the assumed S/N since L2s 
have, on the average, weaker emission lines. For example, as- 
suming S/N=l in the Ka03 sample, I get N(L2)/N(S2)=2.5. 
The samples were selected by their reddening corrected line 
ratios. This introduces only a small difference compared with 
the selection by uncorrected line ratios. 

An additional small sample is taken from Ho et al. 
(1997). These are high quality observations of nearby sources 
observed through small apertures. Here the contribution 
from SB regions is small and can be more easily removed. Us- 
ing the Ho et al. definitions, I choose from the tables only L2s 
and S2s and avoided transition sources (e.g. those marked 
as L2/T2). BH mass estimates for most of these sources can 
be obtained from <r* measurements listed in the literature. 
The archive used for this purpose is LEDA (Paturel et al. 
2003) 



2.2.2 Bolometric luminosity indicators 

Every narrow emission line is a potential indicator. 
However, lines like [N n] A6584 are very sensitive to the gas 
composition (e.g. Groves et al. 2006b) and are not useful 
in this respect. This is not the case for the oxygen and the 
hydrogen Blamer lines that can be used to define four such 
indicators: 

1. The H/9 indicator. In principle, this is the best indi- 
cator because of the very flat dependence of L^H^/Lboi on 
continuum shape and ionization parameter (Fig. [2] and Ta- 
ble [TJ. The normalization of L(H / 3)/Lboi can be obtained 
from type-I observations where the optical AGN continuum 
is directly observed. A major source of uncertainty is the dif- 
ficult deblending of the relatively weak narrow component 
from the strong, broad H/9 line. The fitting procedure takes 
this into account by imposing several criteria, such as assum- 
ing the same profile for the narrow [O ill] A5007 and H/3 lines 
and forcing lower and upper limits on [O iii]/H,3 (see NT07). 
The typical uncertainty in this procedure is estimated to be 
0.2 dex. 

The values of L(H^)/Lboi obtained from the NT07 sam- 
ple are shown in Fig. [3] as a function of [OiiiJ/H^. Also 
shown are means and standard deviations in bins of 0.1 
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Figure 3. Estimated L([0 m])/L bol (black) and L(H,3)/L bc> i 
(red) in 0.1 ^ z $C 0.2 type-I AGN: small points represent individ- 
ual objects from the NT07 sample. Large squares with error-bars 
are means and standard deviations in bins of 0.1 dex in [O Illj/H^. 
The smooth curves are the theoretical models of Fig. [2] (thin lines 
SED-1, thick lines SED-2) where the ionization parameter is re- 
placed by [O Ill]/H^. The lines are shifted down relative to Fig. [2] 
by 1.4 dex. See text for more explanation. 



dex in [O lll]/H,g. The typical standard deviation is 0.25-0.3 
dex and is a combination of the uncertainties in the L(H/3) 
measurements and the scattering in covering factor (see be- 
low). Measurements of L([0 lll])/Lboi in the NT07 sample 
are also shown. As expected from the photoionization cal- 
culations (Fig. [2}, L([0 m])/Lboi increases almost linearly 
with [Om]/H(3 while L(H,3)/Lboi is only weakly dependent 
on [Om]/H,3. 

To make a direct comparison between models and ob- 
servations, I used the calculations shown in Fig. [2] to con- 
vert ionization parameter to [O in] /Tip. For the case shown 
here I chose the calculated L(H/3)/Lboi and L([0 in])/L bo i 
in the SED-1 and SED-2 dusty, solar composition clouds. 
This enables me to plot the theoretical curves of Fig. [5] on 
top of the data in Fig. [3] by shifting the curves down to 
fit the observations. The required shift is 1.4 dex and the 
agreement between model and observations is very good. 
The difference between the vertical scales of Fig. [2] and 
Fig, [3] is a manifestation of the fact that the NLR cover- 
ing fraction is much smaller than unity. The data shown in 
Fig. [3] suggest that for narrow H/3 lines in low ionization 
([Olll]/H,g< 4) type-I AGN, L 6oi ~ 14, 000L(H/3). The ra- 
tio increases smoothly to about 20,000 in high ionization 
parameter sources ([Om]/H/3~ 10). 

The present work uses reddening corrected line intensi- 
ties yet the lines in the NT07 samples are not corrected for 
reddening. The correction factor can be estimated from the 
mean measured Tic/Up in type-II sources. As explained in 
K03, and noted earlier, the amount of reddening depends on 
the stellar population and is higher in smaller D n 4000, large 
SFR galaxies. It is therefore important to split the popula- 
tion into groups, e.g. S2s and L2s. The median H/3 redden- 
ing correction factors for S2s in the Ke06 sample are ~ 3 for 
galactic type reddening and ~ 5.7 for the A -0 ' 7 reddening 
law. The corresponding factors for L2s are ~ 1.6 for galac- 



tic reddening and ~ 2.2 for the A approximation. For 
comparison, the median correction factor in the Bassani et 
al. (1999) sample, for galactic-type reddening, is close to 8. 
The S2 factors are more appropriate since there are very few 
LINERs in the type-I sample I used for the calibration (see 
§3). As explained, narrow emission lines in type-II sources 
suffer, on the average, more reddening compared with type-I 
sources. The additional factor is subject to some uncertainty 
and is estimated to be about 0.1-0.2 dex for H/3. 

Given all these considerations, including the differences 
in extinction between type-I and type-II sources, result in 
the following relation for reddening corrected L(H/3) and the 
A -0 ' 7 extinction law, 



logL boi = log L(Hp)+3.4S+max 



O.,0.31(lo{ 



[Qui] 



0.6) 



For galactic reddening, the factor 3.48 should be replaced by 
3.75. The uncertainty on this ratio depends on the accuracy 
of the scaling shown in Fig. [3] and is estimated to be at least 
as large as the size of the error bars shown in the diagram, 

1. e. about 0.3-0.4 dex. 

The above expression enables me to calculate the mean 
NLR covering fraction. The estimate depends on the amount 
of line photon destruction by internal dust which is already 
included in the photoionization calculations of Fig. [21 As- 
suming this accounts for a factor of 1.5-2 out of the total 
extinction, and a reddening correction factor of about 6 for 
H/3 in S2s, gives a mean NLR covering fraction of 7-15%. 
For galactic type reddening, the covering fraction is smaller 
by a factors of 1.5-2. Obviously, the covering factor changes 
from one object to the next which is probably the main rea- 
son for the large scatter in L(H,3)/Lboi seen in Fig. All 
these estimates are based on the assumption of little or no 
extinction of the optical continuum in type-I AGN which, in 
itself, is subjected to some uncertainty (e.g. Netzer 1990). 

The obvious limitation of the L(H/3) method is the dif- 
ficulty in estimating the SB contribution to the observed 
Balmer lines, especially in large aperture observations. This 
has been discussed, extensively, in K03, K06, Groves et al. 
(2006a), and KH09. K06 define, for each source, a quantity 
which describes its distance from the SB region in the rel- 
evant BPT diagram. Groves et al. (2006a) and KH09 sug- 
gested somewhat different procedures that serve the same 
purpose. The H/3 method is more reliable in cases where 
the SB contributions are small, e.g. extreme S2s with large 
[OtiiJ/H^ or pure L2s. It is also the preferred method for 
small aperture observations, like the Ho et al. (1997) sam- 
ple. 

2. The [Om]A5007 method. Estimates of Lboi based on 
the observed intensity of the [O in] A5007 line are explained 
and discussed in Heckman et al . (2004), Netzer et a. (2006), 
K06, N09, KH09 and several other papers. As demonstrated 
in those references, for high ionization type-I AGN, the scal- 
ing factor between the observed L([Om] A5007) and L boi is 
about 3000-3500. Netzer et al. (2006) show the clear depen- 
dence of this scaling on source luminosity and/or redshift. 
This effect is not taken into account in the present work be- 
cause the type-I and type-II SDSS populations are different 
in terms of the fraction of low ionization (L2) sources and 
I prefer to use, instead, average values for the S2s (see dis- 
cussion below) . Using any [O ill] A5007 luminosity indicator 
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includes, therefore, another uncertainty which is related to 
the source luminosity. 

The use of reddening corrected line intensities reduce 
the conversion factor by a factor of 3-6, depending on the 
extinction law used. For the A -0 ' 7 extinction law I find 
L boi ~ 600L([Om] A5007) in high ionization S2s. This is 
the number used in KH09. For galactic type extinction, 
L bol ~ 1100L([Om] A5007) 

As explained, the major limitation of the [O ill] A5007 
method is the strong dependence of L([0 lll])/Lboi on 
the level of ionization. For a given Lboi, the range in 
L([0 ill] A5007) is more than an an order of magnitude. In 
particular, a constant assumed L([0 m])/Lboi which is cal- 
ibrated from S2 observations severely under-estimate Lboi 
in L2s. This issue is illustrated in the simulations described 
below. 

3. The [O i] A6300 method. The [O i] A6300 line is strong 
enough, in many AGN and SBs, to serve as an additional 
Lboi indicator. The total observed range in [O l]/H Q is about 
1.2 dex, similar to the observed range in [O IIlJ/H^. Thus, the 
expected uncertainty in L([0 l])/Lboi, assuming the Balmer 
lines are the best Lboi indicators, is large. However, the cal- 
culations (Fig. [2} show that the [Oi]/Hq, ratio is not very 
sensitive to the level of ionization of the gas, or its metallic- 
ity, and most of the range in [O l]/H a is likely to be due to 
the changes in SEDs (e.g. Groves 2006a) and/or the column 
density of the NLR gas. There are clearly some trends in 
[Ol]/H a that are not entirely understood. 

Given the theoretical calculations, and the mean ob- 
served line ratios in the Ke06 sample, I estimate Lboi — 
12,000L([Oi] A6300) for S2s and L bot ~ 5000L([O i] A6300) 
for L2s for the A -0,7 reddening law. For galactic extinction, 
the numbers are larger by a factor of ~ 2. 

4. The [O i]/[0 in] method. This is a new indicator which 
is based on the theoretical calculations, the relative weakness 
of [O i] A6300 in SB dominated systems, and simple simula- 
tions that are discussed below. 

The theoretical calculations shown in the right panel 
of Fig. [5] suggest that the [Oi]/[Om] line ratio decreases 
linearly with the ionization parameter. It can thus be used 
to estimate the changes in the level of ionization of the gas. 
This can provide the needed correction for the [O ill] A5007 
method. The calculations also show that the above line ratio 
is not very sensitive to the gas composition, at least within 
the range considered here (1-4 times solar). 

To illustrate this point, I show in Fig. [4] the observed 
[Olll]/H/3 vs. [Oi]/[Oiii] for the high S/N Ka03 and Ke06 
samples. The clear strong trend confirms the similar depen- 
dences of the two line ratios on the level of ionization of 
the gas. As explained, the Ke06 sample was chosen to be re- 
moved, as much as possible, from the SB region. This sample 
is shown by red points and the Ka03 sample by black points. 
I have also added the Ho et al. (1997) S2s (green triangles) 
and L2s (blue triangles). As expected, these small aperture 
observations lie in the Ke06 zone away from the SB region. 

To test the agreement with the theoretical calculations, 
I draw three lines corresponding to three calculated models 
with various assumptions about metallicity, SED and the 
role of radiation pressure (see figure caption). All models 
pass through the Ke06 and the Ho et al. (1997) points ver- 
ifying that the red colored part of the diagram is indeed 
dominated by the AGN contribution. The curves that rep- 




[OI]/[OIII] 



Figure 4. [OiiiJ/H^ vs. [Ol]/[Om] for sources in the Ke06 (red 
points) and Ka03 (black points) samples. The Ho et al. (1997) 
S2s are shown as green triangles and the L2s as blue triangles. 
The black curves represent three theoretical dusty cloud models: 

1. VI SED from Fig.[T]with 3xsolar abundance (thick solid line). 

2. VI SED and constant pressure dominated clouds (thick dashed 
line). 3. SED-1 solar composition constant pressure clouds (thin 
dashed line). 



resent models with different abundances also show that the 
gas metallicity does not affect much the above line ratios. 
The observations and the models suggest the following linear 
relationship in the AGN dominated part of the diagram, 



log 



[QUI] 



= -0.75 log 



_\OI]_ 

[out] 



0.05. 



(2) 



The uncertainties on the two terms, from the fit procedure 
only, is about 10%. Since [Oiii]/H,3 is an ionization parame- 
ter indicator, and since this parameter determines the exact 
L([0 m])/Lboi (Fig. [2}, the two can be combined to derive 
the following expression for the case of a A -0,7 reddening 
law, 

logLboi = 3.53+0.25 log L([0 II J]A5007)+0.75 log L([0 J]A6300) .(3) 

For galactic reddening, the constant 3.53 should be replaced 
by 3.8. The normalization of Lboi is based on eqn.[2]at the 
limit of the largest [OniJ/H^ and the agreement between 
theory and observations in type-I sources (Fig. 01 . The un- 
certainties for individual measurements are due mainly to 
the uncertainty in the slope of eqn.[2j the range in the NLR 
covering factor and the assumed bolometric correction fac- 
tor, BC. The combined uncertainty cannot be smaller than 
- 0.4 dex. 



2.2.3 Simulated BPT diagrams 

To further investigate these issues, I ran simple simulations 
designed to mimic the distributions of SBs, AGN and com- 
posite sources in the [OmJ/H^ vs. [Ol]/H a plane. The sim- 
ulations are shown in Fig. [5] The starting points are the dis- 
tributions shown in K06 Fig. lc and Ho (2008) Fig. 2. These 
are used to define "pure SB" (blue points) and "pure AGN" 
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Figure 5. Simulated BPT diagrams. Pure AGN are marked with 
red points and pure SB with blue points. Left panel: The assumed 
locations of SBs and AGN and observed sources from the Ke06 
sample (small black points). Right panel: Simulated composite 
sources (large black squares) and observed Ke06 sources (small 
black points). The green curve marks the Ke06 AGN-SB division 
line. 

(red points) regions in the left panel of the diagramQ Also 
shown are all sources from the Ke06 sample (small black 
points) and the Ke06 division line separating SB from AGN. 
These simulated SBs and AGN were used to create a sample 
of composite sources that represent different combinations 
of the two in terms of their Lboi (represented by the H/3 
luminosity) and different line ratios. 

I start by randomly mixing the two populations us- 
ing their H/3 luminosity, i.e. for a total luminosity L(H/3), 
a fraction /L(H/3) is contributed by the AGN and a fraction 
(1 — /)L(H/3) by the SB, where / is distributed uniformly 
over the range 0-1.. I then chose, randomly, using the blue 
(SB) and the red (AGN) areas, values for [Oni]/H^ and 
[O l]/H a for the simulated SB and AGN. These are combined 
to form a composite spectrum with composite [O niJ/H^ and 
[Ol]/H line ratios. These simulated sources are shown as 
black squares in the right panel of Fig. [S] The composite 
sources are the only ones analyzed here and the points rep- 
resenting pure AGN and pure SBs in the left diagram are 
only shown to demonstrate their locations. Clearly, there is 
a good agreement between the areas occupied by real and 
simulated sources in the right panel of Fig. [5] 

It is important to note that the density of the simulated 
sources in Fig. [5] is not the same as the density of observed 
sources. Instead, the simulations assume evenly spread val- 
ues of / and evenly spread [O Illj/H^ and [O l]/H Q line ratios 
in the SB and AGN regimes. The diagrams emphasizes the 
problem of identifying the composite sources whose locations 
overlap with the pure SB and pure AGN regions. 

For each of the simulated composite spectra I calcu- 
lated the "real" Lboi of the AGN (i.e. the one obtained from 
/L(H/3)) and compared it with the "estimated" L 00 i, i.e. 

2 Note that B04 assigned a single location to all AGN and KH09 
two representing points, one for S2s and one for L2s. Here I use 
a much large range determined by the observations 
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Figure 6. Estimated Lboi relative to the real (i.e. assumed by 
the simulation) Lboi using the [O III] A5007 (black points) and 
the [Ol]/[Om] (blue points) methods. The upper panels assume 
that the SB contributions to the [O in] A5007 and [O i] A6300 lines 
are perfectly subtracted. The bottom panels make use of the to- 
tal "observed" [O in] A5007 and [O i] A6300 luminosities (i.e. the 
combined AGN-SB line luminosity). The left panels compare the 
two methods for various fraction of the total H/3 flux emitted by 
the AGN. The vertical line marks the 1/3 division line, i.e. the 
SB fraction in left to the line is greater than 2/3. 



the one that would have been obtained by using the var- 
ious Lboi indicators. Four examples are plotted in Fig. Ifj] 
The top two panels compare estimated-over-real Lboi us- 
ing the [O in] A5007 (black points) and the [O l]/[0 ill] (blue 
points) indicators. The two are plotted against the ob- 
served [OlllJ/H^ and the AGN fraction in the composite 
spectrum. This simulation is done under an optimistic as- 
sumption that the SB contribution to each of the lines can 
completely and accurately be removed. As expected, the 
[O ill] A5007 indicator over-estimates the real Lboi in sources 
with large [OmJ/H^ and underestimates it in sources with 
small [Oiii]/H^. The [Oi]/[Oiii] indicator is distributed 
more uniformly, in a narrower band, with a mean of 1.0 and 
a scatter of about 0.15 dex. The bottom panels show similar 
ratios based on the "observed" [O ill] A5007 and [O i] A6300 
line luminosities (i.e. the combined AGN-SB luminosity). 
This results in a larger scatter but the main conclusions are 
unchanged. I have also tested the [O i] A6300 method (not 
shown in the diagram). The results are somewhat inferior 
to the [0i]/[0hi] method but superior to the [Oin]A5007 
method. Obviously, the [O ill] A5007 method looks somewhat 
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Figure 7. Simulated [OinJ/H^ vs. [Ol]/[Oin] diagrams obtained 
from the data shown in Fig. \E\ The left panel shows the entire 
simulated sample and the right panel only sources above the Ke06 
division line. The green points in both diagrams are sources where 
the estimated Lj, ; obtained by the [Ol]/[Om] method deviates 
from the "real" (i.e. assumed in the simulation) Lf, [ by more 
than a factor two. 



Figure 8. A comparison of the various methods for estimating 
L/LEdd in several BH mass and galaxy type (D n 4000) groups for 
Ka03 2 < 0.1 AGN. The histograms show the fraction of sources 
in bins of 0.1 dex in L/LEdd- Black curves are values obtained 
with the [O III] A5007 method, red curves with the H/3 method 
and blue curves with the[0 l]/[0 III] method. The bottom panels 
show L/L^dd distributions in two BH mass groups for all D n 4000 
classes and the top panels two D n 4000 groups with a range of 
-Mb Hi as marked. 



better, and the [O i] A6300 method is more biased, when the 
horizontal axis is replaced by [Ol]/H a . 

A critical issue is the fraction of SB contribution to 
the [O in] A5007 and [O i] A6300 lines that can be tolerated. 
The left panels of Fig. [6] address this issue by showing the 
1/3 division line, i.e. the location where only a third of the 
H/9 line is due to the AGN emission. At this mixing, the 
[O i]/[0 in] method is still producing Ltoi estimates that are 
within a factor 2 of the real values. This number seems to 
be a practical limit when estimating the contamination of 
the [O ill] A5007 line, and hence the uncertainty on Z/ 00 i, in 
optically selected samples. 

Finally I show in Fig. [7] the simulated source distri- 
bution in the [Oiii]/H/3 vs. [Oi]/[Oiii] plane. This plot is 
very similar to the one obtained from the real observations 
(Fig. 2| . There are two parts to the diagram. The left hand 
side shows all composite sources (black points), pure SBs 
(blue) and pure AGN (red). It also shows in green all com- 
posite sources where the [O i]/[0 ill] indicator results in L(, ; 
which deviates by more than a factor of two from the real 
Ltd (12% of the sources). As expected, most of these sources 
lie close to the border line between SBs and AGN. The right 
hand side is a similar diagram for sources in the Ke06 part of 
the diagram. The fraction of green sources here is only 2.5%. 
Thus, using the [O i]/[0 ill] estimator one expects more than 
85% of the sources in the Ka03 AGN region to be within a 
factor of two of the intrinsic Ltoi of the AGN. In the Ke06 re- 
gion, the fraction is more than 95%. Note again that the frac- 
tions quoted depend on the assumed source density across 
the entire plane. Hence, they only apply to samples that are 
distributed with equal density over the assumed AGN and 
SB regions. 



2.2.4 Comparison of different Ltd indicators 

I compared several L/LEdd distributions calculated with the 
various methods. The bottom panels of Fig. [8] show the dis- 
tributions in two BH mass groups, one with 10 7 <M(BH)< 
lO 73 Af and one with 10 s ^M(BH)< 10 8 - 3 M Q . As dis- 
cussed in KH09, lower mass BHs are found, primarily, in 
hosts with younger stellar populations and are, on the aver- 
age, faster accretors. The large BHs are, typically, in hosts 
with older stellar population and lower L/LEdd- The dia- 
gram confirms this finding but shows that the [O ill] A5007 
indicator gives consistently smaller L/LEdd by a factor of 
~ 4. The bottom left panel shows also that the H/3 indica- 
tor over-estimates L/LEdd in lower Mbh sources, because 
of the SB contribution to the line. The bottom right panel 
shows the good agreement of the H/3 and [Oi]/[Oiii] indi- 
cators in high mass BH systems where the SB contribution 
is negligible. 

The upper panels of Fig.[8]compare L/LEdd in different 
groups of the stellar age indicator, D n 4000 (see e.g. K03). 
It suggests that when using the L([0 ill] A5007) indicator, 
L/LEdd is underestimated by a factor ~ 3 in younger sys- 
tems and by a factor of ~ 5 in older systems compared with 
the [Oi]/[Oiii] method. The recent KH09 paper shows a 
powerlaw distribution of L([0 iii])/Mbh in old, large D n 4000 
systems. According to KH09, this is also the distribution 
of L/LEdd in those galaxies. This change in the shape of 
the accretion rate distribution is explained by KH09 as due 
to a different mechanism (stellar mass loss) of mass sup- 
ply to the central BH. The distribution of L/LEdd in the 
1.7<D n 4000< 1.9 group shown here which is based on the 
[Oi]/[Oiii] method does not confirm this idea. Its shape 
cannot be fitted by a powerlaw and its mean L/LEdd is 
larger by a factor of ~ 5 compared with the (translated from 
L([Oiii])/M B h) L/L E dd found by KH09. This is a manifes- 
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Figure 9. Same as Fig.[8]for all Ke06 AGN at all redshifts divided 
into S2s and L2s. 



tation of the failure of the [O ill] A5007 method to estimate, 
properly, L oa i in low ionization AGN. 

Fig. [9] shows a comparison of L/LEdd distributions for 
all S2s and L2s in the Ke06 sample. The two groups, com- 
bined, cover a large range in BH mass, 1O 6 ' 3_9 M0, and have 
different L/LEdd- The diagram shows that the [Oiii] A5007 
indicator is appropriate for the fast accreting S2s (left panel) 
but it underestimates L/LEdd in the slower accreting L2s, 
where the two other indicators agree very well. Again, the 
histograms are not meant to show the real distribution of 
L/Lsdd in all accreting BHs but rather to point to the most 
reliable Ltoi indicators. 



3 DISCUSSION 

Having established [Ol]/[Om] as the best L^oi indicator for 
type-II AGN, I now address two important physical issues: 
the comparison of mass and accretion rates in low redshift 
L2s and S2s and the correlation of L 00 i with SFR and the 
specific SFR (SSFR). Unlike the Heckman et al. (2004) and 
the KH09 papers that address the combined AGN-SB pop- 
ulation, the present work discusses only AGN-dominated 
sources, those where Ltoi is larger than the SF luminos- 
ity. I made no attempt to take into account accreting BHs 
in sources classified as SB galaxies. Such sources fall inside 
the SB region on the Ka03 BPT diagrams and play an im- 
portant role in the KH09 analysis where, in several of the 
diagrams, they outnumber AGN. 



3.1 Black hole masses and accretion rates in low 
redshift S2s and L2s 

A major purpose of this work is to study mass and accre- 
tion rate distributions, and their various correlations, for 
L2s and S2s in the Ka03 and Ke06 samples. Since <r»-based 
estimates of BH mass are known to be reliable in bulge and 
elliptical galaxies, and unreliable in disks and pseudobluges, 
I made an attempt to remove the less reliable estimates. The 
simplest approach is to remove blue galaxies thus avoiding 
disks and pseudobulges (e.g. Drory & Fisher 2007). Since u 



Figure 10. L/L Edi distributions of Ka03 z ^ 0.1 AGN with 
Mbh groups as indicated. The distributions are shown separately 
for red and blue galaxies. 



and r magnitudes are available for the entire sample (SDSS 
"model magnitudes"), f used the u — r method described 
by Baldry et al.,(2004) to separate red and blue galaxies in 
the Ka03 and Ke06 samples. About 75% of the AGN in the 
Ka03 sample and about 85% in the Ke06 sample are found 
to be in red hosts. As expected, the fraction of S2s in blue 
galaxies is larger than their fraction in the entire popula- 
tion. The method is probably too conservative since there 
are bulge galaxies that fall into the "blue" category, espe- 
cially the larger cr* systems. Also, the border line between 
the groups is not as sharp as the expression suggested by 
Baldry et al. All remaining analysis that refers to BH mass 
and L/LEdd in S2s and L2s is based on the red sub-sample 
only. For those cases where I compare type-I and type-II 
properties, I use the entire sample since no red sub-sample 
is available for type-I AGN. 

Fig. [TO] shows L/LEdd distributions for several Mbh 
groups in z ^ 0.1 blue and red hosts in the Ka03 sam- 
ple. The left panel shows the distributions in bins of 0.1 
dex in L/LEdd in the red sub-sample. As seen, the larger 
BHs in the local universe are the slower accretors, similar to 
the behaviour in low and high redshift type-I samples. The 
L/Lsdd distributions are broader than observed in type-I 
sources. This is mostly due to the lack of L2s in type-I sam- 
ples. There is no strong indication for a change in the shape 
of the distribution when going to larger mass BHs. The right 
panel shows the same distributions for the blue sub-sample. 
Here there is no difference between the mass groups reflect- 
ing, most probably, the unreliable Mbh estimates in such 
hosts. Mixing the two sub-samples will cause an increase in 
the intrinsic scatter and will tend to over-populate the larger 
L/LEdd part of the diagram. 

The dependences of L/L E dd on [Oi]/[Oiii] and L bot 
for AGN in red hosts in the Ke06 sub-sample are shown in 
Fig. 1111 The S2s and L2s are marked with different colours 
to demonstrate the smooth transition between the groups. 
Plotted also are L2s from Ho et al. (1997) that probe a much 
lower Ltoi range. The left hand side shows that, without the 
Ho et al. L2s, there is a strong apparent correlation that 
may be interpreted as an indication that the [Oi]/[Oiii] ra- 
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Figure 11. Kc06 AGN in red galaxies divided into S2s (red 
points) and L2s (blue points). The left panel shows L/LEdd vs. 
[O i] / [O III] and the right panel L/L Edd vs - ^bol ■ Large blue points 
denote L2s from the Ho et al. (1997) sample. The two diagonal 
strips in the right panel are two subclasses of the S2s and L2s pop- 
ulations. The top one contain AGN with 7 < log Mbh < 7.3A/0 
and the bottom one AGN with 8 < \ogM BH < 8.3M . 



tio is a good accretion rate indicator. Adding the Ho et al. 
L2s change this conclusion and illustrates the double-nature 
distribution of L/LEdd with respect to [0i]/[0iii]. 

The right hand side of the diagram explores the lumi- 
nosity dependence of L/LEdd- As found in several earlier 
studies, Ltoi and L/Lsdd are strongly correlated over more 
than three orders of magnitudes in the SDSS sample, and 
over more than five orders of magnitude when including the 
Ho et al. L2s. The L/LEdd range in the S2 class is smaller, 
~ 2 orders of magnitudes. The Ho et al. sources fall on the 
continuation of the SDSS correlation. All L/Lsdd and Lf, i 
for L2s would drop by about a factor 5 if the [O ill] A5007 
would have been used instead. Mbh can be considered as 
an additional dimension of the L/Lsdd vs. Ltd correlation. 
This is illustrated by the two parallel bands where the up- 
per stripe show L2s and S2s with Mbh=1O 7_7 ' 3 M0 and the 
lower one Mbh=1O 8_8,3 M0. There are clear and smooth 
transitions between S2s and L2s and between all BH mass 
groups, over more than five orders of magnitude in Lbol and 
no indication for a different sources of mass supply. 

Fig. [T5] shows a comparison of 0.1 < z < 0.2 type-I 
and type-II SDSS AGN. In this case I use the entire Ke06 
type-II sample since the type-I group is not divided accord- 
ing to galaxy colour. The diagrams show L/L B dd vs. L^ i 
(left panel) and L/Lsdd vs. Mbh (right panel). The divi- 
sion between S2s and L2s is the same as in Fig. [11] and the 
type-I sources are shown by small black squares. Ltoi for 
type-Is is measured as explained in N09 (the [O ill] A5007 
method) which is appropriate for such objects. As seen, the 
type-I sources overlap, accurately, the S2s region in both 
diagrams. LINERs are completely missing from the type-I 
sample which results is a much narrower L/LEdd range of 
such sources for all mass BHs. 

The distribution of L/Lsdd in high redshift AGN sam- 
ples has been studied, extensively, in several recent publi- 
cations. Kollmeier et al. (2006), Netzer et al. (2007), Shen 



Figure 12. Left: Same as right panel of Fig. 1111 but for 0.1 < 

z < 0.2 AGN. Small black squares are type-I AGN from NT07. 
Right: the same for L/LEdd vs - A^bh- 

et al. (2007), and Gavignaud et al. (2008), used various se- 
lected samples, with different redshifts and depths, to ad- 
dress this issue. The results of the various studies are quite 
different due to, among other things, the selection methods 
and the way used to estimate Mbh- All those samples, as 
well as other low redshift type-I SDSS samples, clearly do 
not contain LINERs. The reason is the much fainter AGN 
continuum and broad emission lines in such sources which 
makes their detection against the stellar continuum very dif- 
ficult. The L2 fraction in the present low redshift sample is 
as large, or even larger than the S2 fraction yet all the above 
studies completely ignored this population. Obviously, there 
may well be fewer L2s at high redshifts yet it is very unlikely 
that such sources are completely missing. This suggests that 
the overall range of L/L B dd in type-I AGN of all BH mass 
is likely to be much larger than assumed so far. 

3.2 AGN and starburst correlations 

3.2.1 Star formation in type-II AGN 

Most sources in the DR4 sample show clear indications for 
starburst activity and long SF history. This has been dis- 
cussed in K03, K06, Groves et al. (2006a; 2006b), KH09 and 
various other papers. Very detailed discussions, with various 
prescriptions of how to deduce the SFR, are given in B04 
and in Salim et al. (2007; hereafter S07). 

According to B04, the SSFR in galaxies hosting type-II 
AGN can directly be obtained from the SDSS spectroscopy. 
This is achieved by calibrating the SB luminosity in "pure" 
SB systems against the stellar age indicator, D n 4000, us- 
ing emission line luminosities, mostly L(Ha) (see details in 
B04). This can then be applied to all AGN where D n 4000 
has been measured to derive the SFR within the fibre. The 
above methods is subjected to large uncertainties, especially 
in older and redder systems, since the calibration of the 
method is based on D n 4000 measurements in only SB galax- 
ies. In older systems, the typical stellar age can reach several 
Gyrs yet most of the Ha emission is produced in very young 
{few x 10 7 yr) HII regions. This results in a big spread of 
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SSFR for a given D n 4000 (see Fig. 11 in B04). Estimate of 
the errors in using this method are given in B04 Fig. 14. 
For example, at the low SF high stellar mass end, where the 
SSFR is below 10" 11 M©/yr/iW,, the 68 per cent confidence 
level approaches 1.2 dex. The error is significantly reduced 
at higher SFRs. In addition, the SFRs for the entire galaxy 
were deduced by assuming that the SFR within fiber has the 
same dependency on colour as the SFR outside of it. The 
B04 recommendation for estimating SFRs for AGN hosts is 
to use the fibre-based values for both the SSFR and M*. 
These numbers can be found in the DR4 archive and were 
adopted here by using the median-based estimates (see gen- 
eral explanation in the DR4 archive). 

The more recent GALEX-based work of S07 provides 
an independent test of the B04 method. This work uses the 
two GALEX UV bands to provide UV-based SFRs for a 
large number of SDSS galaxies. According to S07, the agree- 
ment between the attenuation corrected (using the Chariot 
and Fall 2000 two component method) Ha-based and UV- 
based SFRs is extremely good for the younger SF systems. 
On the other hand, large deviations can be found in older, 
more massive galaxies, where the B04 method tends to over- 
estimate the SFR. It is however clear from their results (e.g. 
Fig. 19) that some SF activity, at a level of Q.IMq /yr or even 
lower, is observed in many AGN, especially the very slow ac- 
cretors. Inspection of the S07 results (e.g. their Figs. 3, 4 and 
19) suggest that the deviations from the B04 estimates are 
large and the Ha-based estimates are poor for galaxies with 
D n 4000> 1.8. Also, the likelihood of no SF (their "no Ha' 
sources) is much larger for those galaxies with total stellar 
mass exceeding about 10 11,25 Mq. A possible reason for this 
over-estimation of the SFR in large mass, old stellar hosts is 
the identification of a LINER excited Ha line as due to SF 
activity. Given all these uncertainties, I chose to avoid AGN- 
hosts with D n 4000> 1.8 and M* > 10 11,25 when comparing 
SFRs and SSFRs to Lboi and L/Lsdd- 

Having obtained Lsf and SSFR for all AGN (excluding 
those removed by the S07-based criteria), I can now examine 
the correlations of these properties with Lboi and L/LEdd- 
In the following I use the same Mbh groups and the same 
division into S2s and L2s as in the previous sections. The 
analysis is meant to show various correlated properties but it 
is not adequate for estimating global population quantities, 
such as the total mass accretion onto massive BHs in the 
local universe. 

3.2.2 Lsf and Lboi correlations 

Lboi in type-I AGN is readily obtained from optical and 
near-IR spectroscopy. Measuring Lsf in those sources is 
hampered by the strong non-stellar continuum and the in- 
tense emission lines that prevent reliable measurements of 
the narrow Balmer lines and/or D n 4000. The only practical 
way to measure Lsf in luminous, high redshift sources is 
based on infrared (IR) indicators, mostly PAH features in 
the mid-infrared (MIR) and cold dust emission in the far- 
infrared (FIR) . This has been discussed in numerous papers 
(see Sanders and Mirabel 1996 for the earlier work and Lutz 
et al. 2008 for references to more recent publications). 

PAH-based and FIR-based estimates of Lsf in low and 
high redshift, high luminosity type-I AGN are described in 
Schweitzer et al. (2006), Netzer et al. (2007), Schweitzer 



et al. (2008) and Lutz et al. (2008). These papers de- 
scribe Spitzer/IRS spectroscopy of two type-I samples. The 
QUEST sample (Schweitzer et al. 2006) contains 28 low red- 
shift (z ~ 0.1) PG QSOs. About half of the sources show 
clear PAH 7.7 /im emission that is interpreted as a SF sig- 
nature and converted to Lsf using the correlation between 
L(PAH 7.7 (im) and L(FIR). The stacked spectrum of all 
other QUEST QSOs shows the same feature with indication 
of Lsf which is only a factor of ~ 2 smaller for the the same 
AGN luminosity. Lutz et al. (2008) extended this study to 
z — 2 — 3, using Spitzer/IRS observations of strong sub-mm 
QSOs. This provided Lsf for sources of much higher lumi- 
nosity. Lutz et al. (2008) contains also estimates of Lsf for 
additional high redshift AGN observed by Spitzer/IRS. The 
assumption in both studies is that Lsf=vL v at 60/im. As 
shown in Netzer et al. (2007) for the QUEST sample, and 
in Lutz et al. (2008) for the high redshift QSOs, there is a 
strong correlation between Lsf obtained from the FIR and 
L5100 over almost four orders of magnitude in Lboi- 

I combined the present measurements of low redshift 
type-II AGN with the above L(IR)-based estimates for type- 
I sources to show the combined LsF-Lboi diagram for AGN 
dominated sources (Fig. [13} ■ The error bars on Lsf for in- 
dividual sources (yellow lines) are adopted from B04. They 
represent a 34 per cent confidence level. The low luminosity 
type-II AGN fall on the continuation of the same correlation 
as the high luminosity, type-I AGN. The slope of the cor- 
relation is about 0.8 (LsFOcLboi O S ) and LsF~Lboi at about 
Lboi = W 41 ergs s _1 which is outside the range of the sample. 

I note that the line shown here is not a fit to the data 
and there was no attempt to estimate the error on its slope. 
Such a fit could not be reliably obtained not only because 
of the upper limits, but also because the SDSS sources, all 
at the low luminosity end, outnumber the higher luminosity 
sources by a huge factor. It is only meant to show that a line 
of this slope seems to connect all AGN-dominated sources 
with real detections over a large luminosity range. Note also 
that S2s (red points) and L2s (blue points) are well sepa- 
rated but follow the same general correlation. This suggests 
that Lsf follows Lboi in low as well as in high accretion rate 
AGN-dominated sources. 

The correlation in Fig. [TH] may be the result of two se- 
lection effects. One is the omission of SF-dominated systems 
above the straight line and the other the absence of pure 
AGN below the line. The missing sources above the line are 
not due to observational limitations. There are many strong 
SB galaxies in the SDSS sample that show weak or no AGN 
activity and hence are not included in the present sample. 
There are also ULIRGs and SMGs with small or no AGN 
activity. These are also not included in this AGN-dominated 
correlation. The bottom part of high Lboi and very low Lsf 
is more problematic and quite typical of flux limited sam- 
ples. Given the exclusion of AGN hosts with large D n 4000 
and large stellar mass, such sources are not expected to be 
missing at low redshifts since the SDSS type-II AGN sample 
is thought to be complete to about z=0.1. However, the high 
luminosity AGN samples, especially those at the highest lu- 
minosity end, are known to be incomplete and the fraction 
of missing sources with very small SFR is not known. A sce- 
nario which is described below gives at least one possibility 
that does not require such sources. 

The left side of Fig. ll4l shows a suggested time evolution 
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Figure 13. L SF vs. L bo i for AGN. The type-II Ke06 AGN are 
the small points with yellow error bars. The red points are S2s 
and the blue points are L2s. QUEST QSOs from Netzer et al. 
(2007) are shown as large black squares and high redshift QSOs 
from Lutz et al. (2008) as large red squares. Empty symbols with 
lines represent upper limits. The slope of the straight line is 0.8. 



from pure SB (thick line) to a combined powerful SB- AGN 
to to a combined weak SB-AGN. This is a "single event" 
that can occur, in principle, several times in the history of 
a galaxy. The long SF episode starts at time io- Some of the 
cold gas finds its way to the center which results in the onset 
of BH accretion at time ti . The AGN rise time is short and 
lasts until ti . This is followed by an intense AGN phase until 
tz. The diminishing supply of cold gas to the SF regions, and 
the central BH, cause a period where the two fade in parallel 
until ti. The two processes may terminate together or, al- 
ternatively, the AGN fading may last a little longer. In both 
cases, Lsf and Lboi are much smaller at large t compared 
with their peak values. This schematic evolution resembles 
the S07 suggestion of a smooth sequence that begins with 
SF galaxies without AGN and extends at its massive end 
to AGN hosts with different levels of SF. In this scenario, 
weak AGN fall are associated with lower SFR relative to 
the peak activity of both. The strong AGN phase has a tail 
that extends into the domain of quiescent galaxies. AGN in 
early type galaxies that do not show any SF activity are 
not shown. The diagram shows two such scenarios with the 
same SFR. One for a high luminosity AGN (top solid line, 
large L 00 i/Lsf) and one with a weaker AGN (dotted line, 
smaller Lb ;/Lsf). 

The above scenario is translated to an Lsf vs. Lbol 
curve in the right side of the diagram. Here the straight 
line represents the correlation of Fig. 1131 Two pure SF- 
dominated systems with different supply of gas are shown 
as two rising and horizontal lines during times to — t2- The 
fading parts, £3 — ti, representing the decreasing branches 
of both Lsf and Lbol, are shown as lines going down par- 
allel to the main correlation. The regions ti — t^ on both 
curves are obviously more complex. They may involve con- 
stant SF and BH accretion rates (i.e. one point on the curve) 
or periods where both Lsf and Lboi increase along the main 
correlation, perhaps up to a point where L/LEdd^ 1- In this 
scenario, there are no missing objects below the correlation 



Figure 14. A schematic SF-AGN evolution sequence. The left 
part shows the relevant times for the SF (thick line) and AGN 
(thin solid line - luminous AGN, dashed line - low luminosity 
AGN) evolution. The time to — t\ precedes the onset of the AGN 
accretion. The peak (Seyfert) AGN activity occurs between t,2 
and 43 and the long decay, up to tn, is the LINER, low accretion 
rate phase. The right part translates this scenario to an Lsf 
vs. Lboi relationship illustrating two tracks for law luminosity 
(bottom dashed line) and high luminosity (top solid line) events 
(time flags are only marked for the latter). The diagonal line on 
the right represents the observed relationship from Fig. 1131 



line in Fig. 1131 unless the AGN accretion continues into the 
quiescent galaxy part with no SF activity. The left upper 
part of the diagram will be filled by SF galaxies that are not 
included in Fig. 1131 These are found in the SDSS sample and 
in several, high redshift LIG, ULIRG and SMG samples. 

The correlation in Fig. [13] can be translated to a ra- 
tio between the bulge growth rate, g(bulge), assumed to be 
proportional to Lsf, and BH growth rate, g(BH), assumed 
to be proportional to Lboi- For BH radiation conversion ef- 
ficiency of t)bh, and SF radiation conversion efficiency of 
t}sf, the ratio is 



gjbulge) 
9(BH) 
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Vbh/0-1 



tjsf/7 x 10- 



10 42 ergs" 



(4) 



In bulge dominated systems, the time integrated ratio of 
these quantities must equal the local measured value of 
M*/Mbh- This number is larger by a factor of at least 
six compared with the value inferred from eqn. [4] It sug- 
gests that even the slowest accreting BHs in AGN domi- 
nated sources, those at the bottom left part of Fig. 1131 with 
SFR~ 0.1 Mq/yt, are still growing at a relative rate which 
is about six times faster than their integrated cosmic growth 
rate. For the fastest accreting BHs in type-II AGN the ratio 
is larger than 20. The constant (115) in the above equation 
is similar to the one found by Silverman et al. (2009) in their 
analysis of the BH growth rate and SFR in zCOSMOS galax- 
ies (see their Fig. 13). These numbers can be translated, in 
a simplistic way, to the ratio of duty cycles between BH 
and SF activity. It it difficult to assign similar numbers to 
the high redshift population since M*/Mbh is not known at 
high redshift. 

Heckman et al. (2004) studied the M*/Mbh ratio in 
SDSS sources and found good agreement between the to- 



14 Hagai Netzer 



tal accumulated M* and Mbh. Their population average is 
consistent with g (bulge)/ g(BH) ~ 1000. This value should 
not be compared with the one in eq. [4] since the present 
work consider only AGN-dominated sources while Heckman 
et al. (2004) included SF galaxies that outnumber AGN in 
the SDSS sample by a large factor. A good idea of the ratio 
of the two is given in S07 where it is shown that 88% of the 
local SF occurs in SF galaxies and 11% in galaxies that host 
an AGN. This number is consistent with eqn U All these 
factors must be related, in a general sense, to the different 
duty cycles of SF galaxies and AGN. 

3.2.3 Lsf Lboi and stellar ages 

A possible way to test the simple theoretical scenario of 
Fig- 03] is to use the measured D n 4000 as an evolution indi- 
cator. K03 and B04 discussed D n 4000 extensively showing 
that it represents the weighted age of the stellar popula- 
tion. An estimate of this age can be obtained from theoret- 
ical models. For example, the single burst model shown in 
Kauffmann et al. (2003b) gives stellar ages of 4x10 s Mys 
for D n 4000=1.3, 10 9 yrs for D„4000=1.5 and 10 Gys for 
D n 4000=2. Obviously metallicity and various assumptions 
about the nature of the burst will change those numbers. 

Fig. \T5\ shows Lboi and Lsf vs. D n 4000 for various 
groups of AGN. The mean luminosities are calculated in 
intervals of 0.1 in D n 4000 and the error is represented by 
the standard deviation in each bin. While the errors shown 
are relatively small, this is partly due to the choice of data in 
the archive which picks for each D n 4000 the median of the 
SSFR distribution. As shown in B04 Figs. 11 and 14, the 
real distribution is much broader. To include these errors I 
plot in red the 68 per cent confidence intervals on Lsf at 
various D n 4000 as given in B04. As explained, Lsf in hosts 
with D n 4000> 1.8 are highly uncertain and hence are not 
shown. Note also that the diagram is for a given BH mass 
group that contains both S2s and L2s. The S2s are the main 
contributes to both Lboi and Lsf in each D n 4000 group. 

The dependence of both Ltoi and Lsf on the age of 
the stellar population age has been known for some time 
(K03, S07, KH09). However, the great similarity in the shape 
of the two is a new feature. It indicates, yet again, a pro- 
portionality of Lboi and Lsf over a large range of stellar 
populations. This is another manifestation of the LsF-Lboi 
correlation shown in Fig. 1131 The general behaviour resem- 
bles the schematic scenario of Fig. Q3] since D n 4000 is an age 
indicator. The drop at D n 4000^ 1.3 is highly uncertain due 
to the small number of sources in this range. 

Division by the appropriate masses enables a compar- 
ison of L/Lsdd and the SSFR with D n 4000. I plot these 
quantities in Fig. QI)] after converting the SSFRs to the same 
scale as L/LEdd (i.e. multiplying M* by the Eddington lu- 
minosity of a solar mass BH) . The overall shape of the upper 
curve is similar to the previous correlation. It shows a slight 
dependence of L/L B dd on Mbh (the smaller BHs are the 
fastest accretors). The lower curve is only shown for com- 
parison since its shape simply reflect the SSFR-D n 4000 cor- 
relation in the DR4 archive used to obtain the SFRs. The 
short red lines are, again, the B04 68 per cent confidence 
intervals on the SSFR. 

Returning to Fig. 1131 and the various diagrams involv- 
ing D n 4000, one can propose a simple explanation for the 
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Figure 15. Lboi (black symbols) and Lsf (blue symbols) vs. 
D n 4000 for 7.5 < logM BH < 7.8M e z < 0.1 AGN in two sam- 
ples: the Ka03 sample (all AGN, empty squares) and the Ke06 
sample (full squares). Red vertical lines show 68 per cent confi- 
dence intervals on Lsf- 
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Figure 16. L/LEdd and SSFR (normalized to the same units as 
L/L Edd ) vs. D n 4000 for Ka03 z < 0.1 AGN in red galaxies in 
three mass groups. Red vertical lines show 68 per cent confidence 
intervals on the SSFR. 



observed range of Lsf and L/LEdd in a given Lboi bin. The 
first may reflect the combination of early and late type stel- 
lar populations for the same luminosity AGN and the sec- 
ond the combination of BHs with different accretion rates, 
all having the same Lboi- This is a natural consequence of 
the schematic evolution between times t2 and i4 shown in 
Fig. [Ml 



4 CONCLUSIONS 

A critical evaluation of the various ways used to estimate 
Lboi in type-II AGN suggests that some of the previous 
methods, in particular the one based on L([0 in] A5007), 
are problematic. For LINERs, and for AGN in hosts with 
older stellar population and large BHs, this method under- 
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estimates L/LEdd by factors of ~ 5. The paper suggests sev- 
eral alternative methods for estimating Lboi and use them to 
investigate the properties of type-I and type-II low redshift 
AGN from the SDSS sample. The main conclusions are: 

1. The best method for estimating Lboi in type-II 
AGN is based on a combination of L([0 ill] A5007) and 
L([0 1] A6300). The method is appropriate for both high 
(S2s) and low (L2s) ionization AGN. 

2. BH mass and L/Lsdd distributions of S2s resemble those 
of similar luminosity type-I sources. The comparison be- 
tween the two types is based on different ways of measuring 
Mbh and care must be taken in using the <r* method in 
blue galaxies, especially disk dominated and pseudobulge 
galaxies. Most type-II AGN are L2s with very small val- 
ues of L/LEdd and there are very few, if any, such sources 
in type-I samples. This results in a much narrower L/LEdd 
distribution for type-I AGN. 

3. S2s and L2s form a continuous sequence of L/Lsdd with 
no indication for a change in the mode of mass supply. The 
overall range in L/LEdd for the entire population is about 
five orders of magnitude while for S2s it is only three orders 
of magnitude. 

4. There is a clear and strong correlation between Lsf and 
Lid over more than five orders of magnitude in luminosity 
in AGN-dominated systems. The low luminosity, low red- 
shift type-II sources fall on the same correlation as the more 
luminous, high redshift, type-I sources. There is no distinc- 
tion between S2s and L2s all the way down to a SFR of 0.1 
Mg/yr and both follow the same LsF-Lboi correlation. 

5. Lboi, Lsf, L/LEdd and the SSFR all follow, in a similar 
way, the D n 4000 sequence up to D n 4000=1.8, the limit im- 
posed here. Combined with the Lboi-LsF correlation, this 
may give a clue to the sequence of events that leads from 
SF to AGN activity in individual sources. 
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